library(tidyverse)
library(here)
library(glue)
library(sf)
library(fs)
library(mapview)
shp <- here("FKNMS_sample_grid/FlaKeys_Grid.shp")
p <- read_sf(shp)
#p %>% st_set_geometry(NULL)
# dissolve
p_d <- p %>% group_by(subregion) %>% summarize(n = n())
# map
mapview(p_d, zcol = "subregion")
Note that the echo = FALSE parameter was added to the code chunk to prevent printing of the R code that generated the plot.